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We study non-equipartition of energy in granular fluids composed by an arbitrarily large number 
of components. We focus on a simple mean fleld model, based upon a Maxwell collision operator 
kernel, and predict the temperature ratios for the species. Moreover, we perform Direct Monte Carlo 
simulations in order to verify the predictions. 
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I. INTRODUCTION 

The equipartition principle states that the total en- 
ergy of an equilibrium macroscopic system is equally 
distributed among the degrees of freedom of its com- 
ponents. For instance, in a classical gas mixture, equi- 
librium statistical mechanics show that equipartition of 
energy between the components is overwhelmingly more 
probable than a state characterized by different energies, 
i.e. the entropy of the system reaches its maximum at 
the equipartition state [y. This result ensures the oc- 
currence of energy equipartition at the thermodynamic 
limit, and allows to deduce a microscopic expression and 
interpretation for the thermodynamic temperature. In 
this context, temperature is found to be the average ki- 
netic energy of a molecule per degree of freedom. This 
property, which leads to the kinetic definition of tempera- 
ture, reads Te = < K >e, where Nd is the number of 
degrees of freedom in the system, and < K >e is an equi- 
librium average of the kinetic energy. This preliminary 
discussion shows the central role played by equipartition 
of energy in the foundations of equilibrium statistical me- 
chanics, and its deep links with fundamental equilibrium 
concepts such as energy bath or thermal equilibrium. At 
the macroscopic level, equipartition of energy is equiva- 
lent to the zeroth law of thermodynamics, which states 
that if two systems are in thermal equilibrium with a 
third system, they are in thermal equilibrium with each 
other. 

However, the above kinetic definition may also be 
applied in non-equilibrium situations, Tjvb = < 
K >NE- This kinetic temperature, which is now de- 
fined as an average with a non-equilibrium distribution 
in phase space, is a less fundamental quantity than the 
equilibrium temperature, and is equivalent to the average 
agitation of the particles in the system. By construction, 
it does not highlight a priori equilibrium features like the 



equipartition principle. Nonetheless, its inclusion in sys- 
tems where thermodynamics are ill-defined may be rele- 
vant practically, as this quantity still provides meaningful 
information about the variance of the velocity distribu- 
tion and about the local average energy in the system. 
This approach has been applied during the last years to 
describe granular fluids, namely systems composed by a 
large number of agitated grains interacting inelastically. 
In this context, the granular temperature is defined ki- 
netically by: 



T(r,i) 



1 mV'^ 



(1) 



'Electronic address: 
t Electronic address: 



Renaud.Lambiotte@ulg.ac.be| 



lbrenig@ulb . ac. beT 



where V=v — u(r, t) is the random velocity, u(r, t) is the 
local mean velocity and d is the dimension of the system. 
The average is now performed with the one particle 
velocity distribution f{r,v;t). This quantity has been 
shown to give a relevant description of granular fluids, 
for instance by paving the way to an hydrodynamic for- 
malism. In contrast, however, the granular temperature 
exhibits non-equilibrium features, such as an anomalous 
Fourier law for the energy flux PJ , and non-equipartition 
of energy in granular mixtures. This non-equipartition 
phenomenon has been first theoretically predicted in the 
case of binary mixtures by 0] , and successfully verified 
by numerical simulations [j| and experiments |^]. It has 
also been observed in the case of rough inelastic spheres, 
where the translational temperature and the rotational 
temperature are different. It is important to note here 
that energy non-equipartition deserves a careful study 
due to macroscopic consequences such as non-negligible 
corrections to the transport coefficients 0, and its role 
in vertical segregation Q . 

In this paper, we introduce a coherent mean field model 
for low density granular fluids composed by an arbitrary 
large number of components. This model rests on a de- 
tailed analysis for the collision frequencies of the different 
components, and generalizes Maxwell-like models intro- 
duced by Ben-Naim fld| and Marconi ^13] . The values for 
energy non-equipartition are derived analytically in the 
low inelasticity limit, and are successfully compared with 
inelastic hard spheres theory and computer simulations. 
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II. TWO RATES MAXWELL MODEL 

In this paper, we consider a binary system composed 
by smooth inelastic hard spheres in d dimensions, i.e. 
the interactions are instantaneous when the spheres are 
in contact and there is no transfer of angular momen- 
tum. Moreover, we focus on a system that is and re- 
main spatially homogeneous in order to focus on the Ho- 
mogeneous Cooling State (HCS). There are K kinds of 
grains, which are characterized by the following mechan- 
ical quantities: their respective mass mi and diameter ai 
as well as their mutual inelasticity coefficients aij that de- 
scribes the energy dissipation during a collision between 
particles i and j. These different properties discriminate 
the components, which are described at the macroscopic 
level by their concentration Xi = ^ and partial temper- 
ature Ti =< 2 '"2^ >i, where the average is performed 
over species i. By definition, these macroscopic quanti- 
ties are constrained by the relations: ^f^i Xi = 1 and 

^i^i = T. Conservation of impulsion implies the 
following collision rule: 
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where angular integrations have been absorbed into the 
time scale. The integer i goes from 1 to K, and Kij is 
defined by: 

if,,[/(v),g(w)]= / dedw— [/(v').9(w')-.f(v),g(w)] 

(5) 

Because of the Maxwell-like structure of the collision op- 
erators, this set equation leads to a closed system of equa- 
tions for the partial temperatures Ti. Straightforward 
calculations lead to the following expressions: 
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where and are respectively the relative velocities 
Vy=Vi — Vj and positions =ri — of the colliding discs 
i and j . The star velocities are their post-collisional veloc- 
ities and e is the unitary vector along the axis joining the 
centers of the two colliding spheres. In the low density 
limit, by assuming that pre-coUisional correlations may 
be neglected the system is described by a system of 
K coupled Boltzmann equations. In this paper, we use 
mean field methods in order to simplify the mathematical 
structure of the collision operator, by assuming that the 
collision frequency between particles i and j, i.e. propor- 
tional to Vij in the case of hard spheres, is approximated 
by the mean field quantity UiJ: 



where we have introduced the normalized mass ratios 
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We rewrite this set into a more compact 



expression after defining related e^j = I — a.i 



dt 
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The next step consists in introducing the quantities 
Ri = , which measure the departure from energ; 
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equipartition. Obviously, the K quantities Ri, i = 1...K, 
are equal to 1 when equipartition takes place, and i?i is 
always equal to 1 by construction. We rescale the time 

dr — dt-^/^^, thereby considering a time scale propor- 
tional to the average number of collision 1-1 in the sys- 
tem. Let us note that this particular choice is arbitrary. 
This leads to: 




(3) 



Let us stress that this kind of assumption is common 
in order to build Inelastic Maxwell Models @,[l3, [12 ■ 
In the case of binary mixtures, a simplified form has 
been introduced by Ben-Naim and Krapivsky 'lO*! , PiJ = 
-^^Ti{t) + Tj{t). Nonetheless, this approximation does 
not describe correctly the mass dependence of the col- 
lision frequencies and does not give correct relations in 
the case of mass-dispersed mixtures. The assumption (O 
leads to a system of K kinetic equations for the velocity 
distributions fi: 



dt 



Ti 
nii 



V2 



rrii 



rrii 



X [{—2^ij — fiji€ij)Ri + {2fiij — fiij€ij)Rj 



K d-1 

/ nil 



,^^/l + i?,-i,,,(2-.,,) 



x[(-2a*1j - Mjieij) + (2mij - Aiij-eijO^j] (8) 

where one verifies that i?i = 1 as imposed by construc- 
tion. In the sequel, we are interested in the stationary 
solutions for small inelasticity parameters. Therefore, 
we solve ^ by perturbation methods assuming that the 
quantities are all small and of the same order of magni- 
tude. The temperature ratios are written under the form 
Rj = 1 -|- ei?L where e is a small formal perturbation 
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coefficient. This development is based on the fact that 
the zero-th order solution of © is i?^ = 1, Vj, which cor- 
responds to equilibrium equipartition of energy. Given 
these quantities, the stationary first order solution of © 
reads: 
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In a system composed by two components, calculations 
are straightforward and lead to the set of linear equa- 
tions: 
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where i?} = by definition. This expression gives the 
first order deviation to equipartition, in the limit of small 
inelasticity. Remarkably, solution ((Tn|) is exactly equiv- 
alent to the solution obtained by Garzo and Dufty 
from the inelastic hard spheres Boltzmann equation in 
the low inelasticity limit. This equivalence, together with 
the mathematical simplicity induced by the mean field 
model, seem to indicate that the TRMM is a relevant 
model for granular mixtures, which could be used, for 
instance, in order to probe the applicability of hydrody- 
namics in these systems, or to derive the temperature 
ratios in mixtures composed by an arbitrary number of 
components, as shown below. 



III. ARBITRARY NUMBER OF COMPONENTS 

In order to treat systems composed by a large number 
of species, it is helpful to rewrite Q into the canonical 
form: 



K 

AiRl+ B,jR]^C, 
where the coefficients read explicitly: 



(11) 
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(12) 



Consequently, the whole dynamical problem is reduced 
to the inversion of the matrix M defined by: 



(13) 



Unfortunately, this problem is practically not trivial 
when the system is composed by a large number K of 
components (inversion of a K-1 matrix) and leads to an- 
alytical but intractable expressions for the Rj. In the 
following, we prefer to focus on a particular case which 
enormously simplifies the inversion of and whose so- 
lutions Rj have a simple analytical form. Namely, we 
consider a system where all particles have the same mass 
nij = m and the same diameter uij = cr, but whose 
inelasticities vary. Remarkably, in that case, all the 
non-diagonal terms Bij vanish exactly Bij = 0, so that 
the general solution of (|12|l reads: 



R] = 



K 



/■(eij 



(14) 



This solution is instructive for several reasons. First, 
let us insist on the fact that it is the first expression 
for non-equipartition of energy in multi-component mix- 
tures. This solution clearly shows that the system be- 
haves qualitatively in the same way as in a binary in- 
elastic mixture, namely the partial temperatures Ti do 
not reach asymptotically the same value, but that they 
remain proportional in the long time limit. This be- 
haviour justifies the derivation of granular hydrodynam- 
ics for complex mixtures. Another nice feature of H14|) 
comes from its very simple expression, which allows to 
recover the solution p(J|l when K = 2, and which is easily 
generalized when the system is composed by an infinite 
number of species, namely when there is a continuum of 
inelasticities in the system. In that case, instead of the 
concentrations Xi used in order to characterize the com- 
ponents, we use the concentration density p(t), which 
represents the density of species l. The probability nor- 
malization reads: 



(15) 



dip{L) = 1 
Relation H14|) generalizes into: 



di- P(t )(e„' - e^,') 



(16) 



where r is the arbitrary chosen reference species to which 
the temperatures are compared, R{l) = 

We verify the validity of (|14|l by performing two kinds 
of computer simulations for the multi-component mix- 
ture. The latter were performed for particles whose 
masses and diameters are equal, m^ = 1 (7^ = 1 in order 
to compare the results with (|14|) . First, we make DSMC 
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FIG. 1: Asymptotic temperature ratio, as a function of the 
inelasticity of each species, in a system where amin ~ 0.8 
(a) and amin ~ 0.4 (b). The solid line corresponds to the 
theoretical predictions USt) . while the data points correspond 
to DSMC simulations of TRMM and IHS. 



which is a reasonable assumption for the inelasticity of 
collisions between particles i and j, and which allows to 
characterize each species i by their sole inelasticity an , in- 
stead of the vectorial quantity {an, ...^au, ...,a]\[i). This 
assumption, which also implies that = ^^"j^i simpli- 
fies into: 

1 ^ 1 

= ^(ei-e.) (18) 

In the simulation results presented below, we define a 
minimum and a maximum inelasticity in the system, 
amin and amax respectively. By definition, we chose 
ail = amin- Then, we fill uniformly the interval 
[amin,amax[ , 1-6. a^i = a^in -I" (i- l)(5a, where the quan- 
tity Sa is defined by Sa — . Let us stress that 
this uniform distribution in the interval [amin , amax [ is 
an arbitrary choice, and that more general compositions 
can be considered without any analytical nor computa- 
tional additional difficulty. In this work, we also perform 
DSMC simulations of the true set of inelastic Boltzmann 
equations for the same mixture, where no approximation 
is made to simplify the collision operator, thereby testing 
the validity of the TRMM. In figure ^ we present results 
in the small inelastic limit amin = 0-8, that show an ex- 
cellent agreement with (|18f) . and in the high inelasticity 
limit amin = 0.4, for which only small discrepancies from 
the predictions may occur, i.e. deviations do not exceed 
15%. 



IV. CONCLUSION 



simulations of the set of kinetic equations Q , namely we 
apply the standard DSMC algorithm, where the proba- 
bility of collisions of a pair is taken to be proportional 
to ^JTT+Tj- In the simulations, we considered system 
composed by a large number if = 50 of components. 
Moreover, each species is composed by the same number 
of particles, Ni = Nj. This implies that the concentra- 
tions are equal to Xi = The species are discriminated 
by their inelasticities an , which are all assumed to be dif- 
ferent, and we chose the reference component to be the 
one with the lower inelasticity: an < an, i > 1. This 
arbitrary choice suggests that the quantities Ri should be 
> 1 in the long time limit, given the fact that other com- 
ponents dissipate less energy than the reference compo- 
nent does. One should note, however, that the choice for 
the cross inelasticities can still be arbitrarily chosen 
provided they respect the symmetry relation = a^j. 
In the following, we choose the following expression for 
the cross inelasticities: 



In this paper, we have focused on the non-equipartition 
of energy in inelastic gases composed by a large number of 
species. We have used mean field approximations, in or- 
der to simplify the set of Boltzmann equations, and to fo- 
cus on an analytically tractable problem. By doing so, we 
have defined the Two Rates Maxwell Model, that we have 
solved formally in the small inelasticity limit, thereby de- 
riving the explicit values for energy non-equipartition in 
systems where components have the same mass and di- 
ameter. These predictions have been verified by simula- 
tions for systems composed by 50 components, that show 
that the TRMM is a relevant model, even in very inelas- 
tic systems. At this point, it is important to note that 
the influence of the shape of the velocity distributions 
has not been taken into account by the mean field mod- 
eling. This is due to the Maxwell-like kernel of the colli- 
sion operators, that leads to a closed set of equations for 
the partial temperatures, thereby neglecting influence of 
higher velocity moments [l3|- Nonetheless, this influence 
is usually weak and has no quantitative effect when the 
system is weakly inelastic. That was shown in the small 
inelasticity limit, both by the DSMC simulation results 
and by analytical comparisons with the exact inelastic 
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Boltzmann equation. As a consequence, the TRMM is 
an ideal candidate in order to apprehend more complex 
phenomena in granular mixtures, such as hydrodynamics 
for a large number of components, or the detailed study 
of systems composed by a continuum species. 
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